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Abstract 



The coupled cluster method has been applied to the eigenvalue problem lattice Hamil- 
tonian QCD (without quarks) for SU(2) gauge fields in two space dimensions. 

Using a recently presented new formulation and the truncation prescription of Guo et 
al. we were able to compute the ground state and the lowest 0"*"-glueball mass up to the 
sixth order of the coupled cluster expansion. 

The results show evidence for a "scaling window" (i.e. good convergence and Constance 
of dimensionless quantities) around (3 = A/ ^ ?>. 

A comparison of our results to those of other methods is presented. 

1 Introduction 

It is the purpose of this paper to present our first results of an attempt to compute the spec- 
trum of a lattice gauge field theory within the Hamiltonian formulation. Our computational 
framework is the coupled cluster method which has been presented and discussed in detail in 



We have applied our method to the lattice version of the 2+1 dimensional SU (2) Yang-Mills 
theory which is the simplest non-trivial non-abelian lattice gauge field theory. We consider this 
model as an important test case for controlling our computational scheme. Although there is no 
experiment to compare with, there exist for this model very reliable results for the ground state 
energy [0] and for glueball masses which allow a critical test of the coupled cluster predictions. 

2 Computational Scheme 

Within the Hamiltonian formulation, the computational framework is given by the Kogut 
Susskind theory which has been discussed in detail in Ref. []T|. In summary, the structure is 
the following: 

The Kogut Susskind wave functions ^{U) = "^{Ui, ..,Un) depend on link variables Ui (/ = 
1, ..,N) which are elements of the gauge group SU{2). N is the number of oriented links in 
a 2-dimensional lattice and is related to the chosen finite volume. Our many-body techniques 
allow us to perform the final calculations in the infinite volume limit N ^ oo. 

The idea of the coupled cluster method is to reformulate the eigenvalue problem = E"^ 
as equations for the ground state correlation function S and for the excitation operators F, 
following from the ansatzes \l/o(f^) = exp S{U) for the ground state and "^{U) = F{U) exp S{U) 
for excited states. 

lfH= l/{2a){J2ik9^Efk-'^/9^y) is the Kogut Susskind Hamiltonian (we use the notations 
{Ik) for link-colour quantum numbers, {k = 1, ..3), Ei^ for the momentum operators conjugate 
to Ui and V for the plaquette term.) this yields the non-linear equation 
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for S and the linear equation 



J2{FiMk + 2SikFik) = '^{E-Eo)F (2) 
Ik 9 

for the excitation operators F. fik abbreviates "hnk variable derivatives" of functions f{U): 

fik = [Elk, f] , fikik = [Elk, [Elk-, /]] • (3) 

This (rigorous) form of the eigenvalue problem guarantees manifestly the correct volume 
dependencies of the ground state energy Eq and the excitation energies E — Eq (see Ref.|l|]). 

Within this paper we will restrict ourselves to excitations corresponding to glueballs defining 
a trivial representation of the lattice Euclidean group (O"*" glueballs with momentum zero). In 
this case we may write 

F = HoFi^t S = Ho^int (4) 

where Ho is the projection operator on states with vanishing (lattice) momentum, angular 
momentum and parity plus. 

A calculation for general representations of the lattice Euclidean group is possible (see||l|) 
but technically more involved. As a first test for our method we have chosen the simplest case. 
Our results encourage us to try the generalization and we are doing this now, but we will defer 
the presentation of results to a forthcoming publication. 

The essential point of the coupled cluster method is to expand the intrinsic functions ^int 
and Fint with respect to a gauge invariant, linked, standardized basis Xa in the form 

SintiU) = S^x'^iU) FUU) = E F^x'^iU) (5) 

a a 

and to define approximations by truncations of this basis. 

Introducing the constant function via = 1, the "plaquette function" by putting IIoX^ = 
AV and using the strong coupling structure Y.iaXfaia — ^aX°'i the coupled cluster equations (|l|) 
and (|^) become equations for the coefficients Sa and F^ 

(^aSa + YCl{S)Sp = TT^^al + Y^^SaO (6) 
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e^F^ + 2Y,C^{S)Ff, = '^(E-Eo)F^ (7) 

with cf(5) = . 
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Here, the coupled cluster matrix elements C^^ are given by 



= l{ea-ep-e,)j:cli (8) 



u 



where the numbers c^l related to the action T{u) {u E Ge) of the lattice Euclidean group 
G^; on the basis x° and it's products by 

E x'T{u)x'^ = J:c^IT{u)x'' (9) 



Hereby, only those cases have to be considered where the functions ^^nd T{u)x^' have a 
common hnk variable. Details are given in 

The basis which we use for our calculations is systematically generated from and 
given by the following prescriptions. 

Introduce first a set of "generic" functions A^*^ given by linked, standardized 5-fold plaquette 
products. They have the structure 

kf = x'nu2{5,k))x\..T{us{5,k))x' ; k = l,..,ns (10) 

and we define no = ni = 1, A^"*^ = A^"*^ = x^- 
5{= 0, 1, 2...) is the "order" of the function A^'^. 

For our two-dimensional SU(2) case we have ns = 1,1,2,4,12,35,129 up to sixth order, 
respectively. The functions (0) are characterized by simple loop patterns exemplified in |l|. 

The essential property of this set of functions is that from them one can construct an 
orthogonal (in the limit 6 —>■ oo complete) basis by acting with suitable Casimir operators 
of the lattice orthogonal group on A^'^ and diagonalizing the corresponding Casimir matrices. 
This basis is then (iteratively) "ordered" by 6 because the Casimir operators do not increase 6, 
i.e. each element x" has some order 6{a). 

Up to order six, one obtains in this manner subspaces of the function space of link variables of 
dimensions (1, 2, 5, 15, 84, 557, 4972), respectively. The relation to the simple generic functions 
(p!0|) , which is set up by diagonalizing the Casimir matrices numerically, allows then to compute 
also the coefficients c^!^ (see |jl|). 

When constructing the basis Xa in the way described, one has to take care of possible 
linear dependencies between the generated functions. In previous investigations [^, ^, |^, |l| 
this problem was solved by exploiting the Cayley Hamilton relation between matrices. We 
have used in this connection a much simpler procedure: If a set fi{U), .., fn{U) contains only 
m {m < n) independent functions, the matrix fi{U''){i,k = l,..,n) has for suitable fixed 
variables ([/^, .., f/") exactly the rank m. Our experience is that statistically chosen variable 
sets (f/^,. .,[/") G SU{2) are suitable in this sense. The linear relation between the functions 
fi is then easily constructed and the dependent functions can be eliminated. 

Having determined in this way all ingredients of the coupled cluster equations (||) and (|^), 
we still have to define a truncation prescription. Up to now, we have used in this connection 
the proposal of Guo et alp which yields actually the simplest set of equations. In this case 
one puts in the order 6 

C;"' = for 5(ai) +5(^2) > 5 (11) 

3 Results and Discussion 
3.1 The Ground State 

It is a special feature of the Hamiltonian formulation that is provides the energy and the wave 
function of the vacuum state which is the ground state of the Hamiltonian. 

Standard lattice Monte Carlo calculations do not give results for the vacuum energy density. 
There exist, however, computations within the strong coupling expansion |Q and very reliable 
Green's function Monte Carlo results 0. 

In Figure 1 we compare our coupled cluster results up to sixth order to the results of the 
other methods. As coupling variable we use the standard expression j3 = 'i/g^. We see that 
we have a good quality and convergence up to /? ~ 4, but our method breaks down for large {3 
where the Green's function Monte Carlo is still valid. 



An important feature of our results is that the vahdity of the coupled cluster method clearly 
goes beyond the range of the strong coupling expansion which breaks down at /3 ~ 3. (Our figure 
gives the 18th order of strong coupling perturbation theory!) We consider this improvement as 
a necessary condition for obtaining continuum limit physics. 

We should stress that the determination of the vacuum operator S{U) - by solving (|^) 
iteratively - turns out be be much simpler and faster than that of the excitation operator F{U). 
In this sense we do not see any special difficulty with the fact that the Hamiltonian formulation 
also involves the determination of the vacuum state. The coupled cluster formulation deals 
with this problem apparently quite effectively. 

Note also that assuming a trivial (strong coupling) vacuum (i.e. S{U) = 0) makes no sense 
within our framework, this would only yield the completely unphysical strong coupling limit 
spectrum. 

In this sense our framework seems to be quite orthogonal to the light front formulation 
which is based on the assumption that the simplicity of the vacuum could be helpful for the 
determination of the spectrum. 

3.2 The Glueball 

The solution of (|^) with the lowest eigenvalue has the interpretion of an (approximate) O"*" 
glueball mass. 

Within the standard lattice Monte Carlo method, a reliable determination of this glueball 
mass has been achieved by TeperQ. 

A comparison to our Hamiltonian results involves a rescaling of the standard "Euclidean" 
coupling qe relative to the Hamiltonian coupling g given by 

/3e = /3 + .077 + O(/) (12) 

and a rescaling of the Euclidean masses Me relative to the Hamiltonian masses M ( "velocity 
of light correction" ) 

Me = il + 0.08Ag'^ + 0{g^))M (13) 

(See Ref.g]). 

According to Teper's results, the region of asymptotic scaling starts at ~ 3 for the string 
tension and for a somewhat larger value of /3e for the glueball masses. 

Since our ground state results are reliable up to /3 ~ 4, one may hope for a scaling window 
for our glueball predictions in this range. 

A direct comparison to Teper's results is difficult in this /3-range because the lowest order 
relations (12) and (13) are not sufficient for this case. These corrections are large in the strong 
coupling regime itself because the strong coupling expansions are completely different in the 
Euclidean and Hamiltonian formalisms; the Euclidean expansion has a logarithmic singularity 
for (3 = which is not present in the Hamiltonian framework. 

Therefore, the best test for a window displaying continuum physics is the computation of 
dimensionless quantities which should become independent of f3. Since we have (up to now) 
only one glueball at our disposal (higher 0+-glueball masses turn out be be rather unreliable), 
we have determined within the same coupled cluster approach the "radius" R of the glueball 
defined by 

R' = (E iFal'lx'-l'Nixn) /(E \Fa\'\xc^\') (14) 

\ a /a 

where A^(x") is the number of different plaquettes occuring in the corresponding generating 
function ([I0|). (The Haar measure norms were calculated using a Monte Carlo method.) 



The result for dimensionless quantity RM is given in Figure 2. There is apparently evidence 
for a scaling window around /5 ~ 3. 

In Figure 3, we also give our results for the absolute values of the glueball masses. Most 
interesting is a comparison to the predictions of the (Hamiltonian) strong coupling expansion [0] 
and of corresponding extrapolations (ELCE-method)[0 which are based on the same Kogut 
Susskind Hamiltonian. 

While the strong coupling expansion breaks down at (3 ~ 2.1 (seen in Figure 3 from the 
28th order of this expansion) our results show a clear convergence (up to order six) up to /5 ~ 3. 
Thus we are able to make reliable calculations in the range which is beyond, but not too far 
beyond the radius of convergence of the strong coupling expansion. This is actually the typical 
range of the reached scaling window in most Monte Carlo lattice QCD calculations. 

Comparing to the ELCE method, we see a clear discrepancy already in the range /5 ~ 3. The 
ELCE extrapolation seems to overestimate the mass values and it is difficult to understand how 
this result will converge to the asymptotic glueball mass, given by Teper as M{g 0)/2g^ = 3.2 
with a 1% error. 

The coupled cluster method, on the other hand, appears to give a prediction that does 
not contradict to this continuum limit value, although a direct comparison to Teper's result is 
hampered by the lack of good knowledge of the rescaling corrections between Hamiltonian and 
Lagrangian formulation for (3^3. 

Summarizing we conclude that the coupled cluster method appears to be able to produce 
predictions within Hamiltonian lattice field theories in a range of the couplings where asymptotic 
scaling is already valid. Hereby, the simultaneous computation of the vacuum state poses no 
special difficulty, it may be reliably determined in the region where observables show a scaling 
behaviour (and beyond, i.e. up to /5 ~ 4 in our case). Also the numerical effort is much less 
for the vacuum than for the diagonalization problem of the glueball masses. 

Clearly, a more decisive test of our computational framework will be our predictions for 
glueball masses with more general (lattice) angular momentum quantum numbers and the 
comparison of the emerging mass ratios to those of Teper. Calculations in this direction are 
now being undertaken. 

We would like to thank Bernhard Metsch for many valuable discussions. D. S. also acknowl- 
edges a useful exchange of ideas with Chris Hamer, Helmut Kroger and Xiang Luo. 
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Figure Captions 



Figure 1: Coupled cluster predictions of this work for the ground state energy density in 3rd, 
4th, 5th and 6th order, compared with the 18th order of the strong coupling expansion and 
the results of the Green's function Monte Carlo method [0. The upper part magnifies the low 
(3 range of the results to show the quality of the convergence. 



Figure 2: The dimensionless product of glueball mass and radius MR in 4th, 5th and 6th order. 



Figure 3: Mass M of the 0^ glueball in 4th ( ), 5th ( — ) and 6th ( ) order in comparison 

with the 28th order of the strong coupling expansion (lower dot-dashed line) and an extrapo- 
lation using the ELCE method |^ (diamonds with error bars). The upper dot-dashed line gives 
the Monte Carlo results of Teper 0] using a lowest order Euclidean-Hamiltonian rescaling given 
by eqs. (12) and (13). The Monte Carlo errors (1-3%) are small compared to the uncertainty 
of these rescaling corrections (see text) and are therefore not indicated in the figure. The insert 
shows a magnification of the crucial part between (3 = 2 and /5 = 3. 
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